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Severe acute respiratory syndrome coronavirus (SARS-CoV) infection can cause the development of severe 
end-stage lung disease characterized by acute respiratory distress syndrome (ARDS) and pulmonary fibrosis. 
The mechanisms by which pulmonary lesions and fibrosis are generated during SARS-CoV infection are not 
known. Using high-throughput mRNA profiling, we examined the transcriptional response of wild-type (WT), 
type I interferon receptor knockout (IFNAR1 - - ), and STAT1 knockout (STAT1 - - ) mice infected with a 
recombinant mouse-adapted SARS-CoV (rMA15) to better understand the contribution of specific gene 
expression changes to disease progression. Despite a deletion of the type I interferon receptor, strong expres¬ 
sion of interferon-stimulated genes was observed in the lungs of IFNAR1 -7- mice, contributing to clearance of 
the virus. In contrast, STATl _/_ mice exhibited a defect in the expression of interferon-stimulated genes and 
were unable to clear the infection, resulting in a lethal outcome. STAT1 -7- mice exhibited dysregulation of 
T-cell and macrophage differentiation, leading to a T H 2-biased immune response and the development of 
alternatively activated macrophages that mediate a profibrotic environment within the lung. We propose that 
a combination of impaired viral clearance and T-cell/macrophage dysregulation causes the formation of 
prefibrotic lesions in the lungs of rMA15-infected STAT1 -7- mice. 


The severe acute respiratory syndrome coronavirus (SARS- 
CoV) is a highly pathogenic respiratory virus that emerged in 
the human population to cause a global epidemic in 2002 and 
2003. The virus infected roughly 8,000 people and resulted in 
approximately 800 deaths, with high case mortality in popula¬ 
tions of the elderly (9). In humans, SARS-CoV infection 
caused mortality due to development of acute respiratory dis¬ 
tress syndrome (ARDS) (19, 23). The lung pathology in SARS- 
CoV-infected patients is believed to be caused by host immune 
dysregulation as well as by viral factors (1, 28, 63). 

ARDS, with a mortality rate of approximately 50%, is the 
most severe form of acute lung injury and affects more than 
one million people per year worldwide (32, 77). No effective 
therapies for ARDS currently exist, making it a major chal¬ 
lenge for critical-care medicine. The first phase of ARDS is the 
exudative phase, presenting with diffuse alveolar damage 
(DAD), inflammation, proteinaceous edema, hyaline mem¬ 
brane formation, and severe hypoxia. This may progress to a 
second phase of ARDS, the organizing phase, characterized by 
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pulmonary fibrosis. In addition to SARS-CoV, viruses that can 
cause ARDS include influenza virus, respiratory syncytial vi¬ 
rus, adenovirus, varicella zoster virus, human metapneumovi¬ 
rus, and hantavirus (5, 12, 34, 37, 38, 47, 69). Mice infected 
with select strains of SARS-CoV or influenza virus serve as 
models for the initial exudative phase of ARDS (63, 64, 80), 
but these animals do not develop pathology consistent with 
later ARDS phases, and there are few infection systems that 
serve as models of pulmonary fibrosis (75). Instead, most stud¬ 
ies focused on pulmonary fibrosis utilize the bleomycin model, 
which does not recapitulate the effects of a viral inducer of the 
disease (50). 

Several studies have shown the importance of interferon- 
and STATl-dependent signaling for reducing viral replication 
as well as for limiting pathogenesis in animal models. Infection 
of type I interferon receptor knockout (IFNAR1 -7- ) animals 
with highly pathogenic H5N1 influenza virus results in accel¬ 
erated mortality compared with wild-type (WT) animal results 
(71). Mouse hepatitis virus, another coronavirus, causes lethal 
disease in type I interferon receptor knockout mice but not in 
wild-type mice (7). Likewise, STAT1 7 “ mice exhibit dramatic 
defects in the immune response, leading to mortality in re¬ 
sponse to Listeria monocytogenes , vesicular stomatitis virus, 
and herpes simplex virus infections (46, 58). In vivo treatment 
with interferon reduces SARS-CoV replication and lung pa¬ 
thology in macaques, suggesting that the innate immune re¬ 
sponse does protect animals from SARS-CoV pathogenesis 
(29). 
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Previous gene expression studies of SARS-CoV in model 
systems have provided some insight into the mechanism of 
pathogenesis. Initial studies using macaques infected with sub- 
lethal inoculums revealed a strong innate immune response 
that corresponded with peak viral replication and a later pro¬ 
liferative response associated with increased expression of cell- 
division genes and the healing of the damaged tissue (17). 
Further studies of lethal SARS-CoV infection in aged mice 
revealed an accelerated upregulation of acute-phase response 
genes associated with ARDS, including the cytokines interleu¬ 
kin-1 beta (IL-1(B), IL-6, and tumor necrosis factor alpha 
(TNF-a) (63). However, infection with clinical isolates of 
SARS-CoV did not cause significant lung disease in young 
mice and resulted in delayed upregulation of immune response 
genes compared with the results seen with aged mice. To 
enable routine studies of SARS-CoV pathogenesis, mouse- 
adapted SARS-CoV strains have been developed using serial 
passage techniques (16, 52, 62). The first of these strains, 
designated rMA15 (representing a recombinant mouse- 
adapted SARS-CoV), contains 6 mutations with respect to the 
Urbani clinical isolate of SARS-CoV (62). Both young and 
aged mice infected with rMA15 develop lung pathology, and 
several strains of mice succumb to infection. Using this virus, 
we have been able to infect young mice and elicit ARDS 
pathophysiology (61). 

Our previous work (24) revealed that young type I interferon 
receptor knockout (IFNAR1 /_ ) mice did not exhibit in¬ 
creased pathogenesis as a result of rMA15 infection compared 
to wild-type (WT) mice. However, rMA15 infection of young 
STAT1 -- mice resulted in increased pathology, including ev¬ 
idence of pulmonary fibrosis-like lesions and mortality by 9 
days postinfection (dpi). In this work, we extended the initial 
biological characterization of rMA15 infection in WT, 
IFNAR1 /_ , and STAT1 /_ mice by using microarray profil¬ 
ing, immunohistochemistry, and cell-sorting techniques to ex¬ 
amine the molecular mechanisms governing severe end-stage 
lung disease. We found that STAT1 /_ mice infected with 
rMA15 exhibited a shift to a T H 2 response, with expansion and 
activation of alternatively activated M2 macrophages, resulting 
in fibrosis-like protein deposition in the lungs. 

MATERIALS AND METHODS 

Mouse infection. 129S6/SvEv (abbreviated here as “129”) wild-type and 
STAT1“~ mice (catalog number 002045-M-F) were obtained from Taconic 
Farms (Germantown, NY). IFNAR1 mice bred on a 129SvEv background 
were a gift from Mark Heise (University of North Carolina [UNC]-Chapel Hill) 
and were bred at the University of North Carolina mouse facility (Chapel Hill, 
North Carolina). All animal studies were conducted using animal biosafety level 
3 laboratories and Sealsafe cages with HEPA filters, and personnel wore per¬ 
sonal protective equipment, including Tyvek suits and hoods as well as positive- 
pressure air respirators with HEPA filters. Ten-week-old mice were anesthetized 
with a ketamine (1.3 mg/mouse)-xylazine (0.38 mg/mouse) mixture administered 
intraperitoneally in a 50-pl volume. Mice were infected intranasally with either 
phosphate-buffered saline (PBS) alone (Invitrogen, Carlsbad, CA) or 1 X 10 5 
PFU rMA15-PBS. A total of 5 infected mice from each strain were euthanized 
at days 2, 5, and 9 postinfection (dpi) for virological, histological, and microarray 
analyses. Mock-infected animals from each of the strains were euthanized at 5 
dpi. Samples from the right lung were reserved for viral titers and histology, and 
the entire left lung was reserved for microarray analyses. Animal housing, care, 
and experimental protocols were in accordance with all UNC-Chapel Hill Insti¬ 
tutional Animal Care and Use Committee guidelines. 

RNA preparation and oligonucleotide microarray processing. The entire left 
lung of each mouse was dissected and homogenized using Trizol and a Magnal- 


yser system (Roche) according to manufacturer instructions. RNA was further 
purified using RNeasy columns (Qiagen) according to manufacturer instructions. 
RNA samples were spectroscopically verified for purity, and the quality of the 
intact RNA was assessed using an Agilent 2100 Bioanalyzer. For each set of 
treatment conditions, 3 of the 5 RNA samples exhibiting the highest RNA 
integrity number (RIN) determined using the Bioanalyzer were used for mi¬ 
croarray analysis. cRNA probes were generated from each sample by the use of 
an Agilent one-color Quick-Amp labeling kit. Individual cRNA samples were 
hybridized to Agilent mouse whole-genome oligonucleotide microarrays (4 X 44) 
according to manufacturer instructions. Samples from individual animals were 
not pooled to enable examination of animal-to-animal variations as part of the 
data analysis. Select samples were hybridized a second time (n — 2 technical 
replicates) to verify the quality of the process. Slides were scanned with an 
Agilent DNA microarray scanner, and the resulting images were analyzed using 
Agilent Feature Extractor software. Data were warehoused in a Katze Labkey 
system (Labkey, Inc., Seattle, WA) and analyzed using Resolver 7.1 software 
(Rosetta Biosoftware, Seattle, WA), and Spotfire DecisionSite for Functional 
Genomics 9.1 software (Tibco Spotfire, Somerville, MA). Primary data are avail¬ 
able (http://viromics.washington.edu) in accordance with proposed Minimum 
Information About a Microarray Experiment (MIAME) standards. 

Analysis of microarray data. An initial analysis of gene expression in each 
mouse line was performed as a function of time (kinetic analysis). Genes were 
selected based on textbook analysis of variance (ANOVA) of microarray inten¬ 
sity scores, with a significance criterion of P < 0.05. The following infection 
groups (n = 3 per group) were compared for WT, IFNAR1 -7- , and STATl -7 ' 
mice, for a total of 9 gene lists: mock and 2 dpi, 2 dpi and 5 dpi, and 5 dpi and 
9 dpi. These gene lists were reduced by requiring a change of at least 1.5-fold 
between the group averages. Genes that met these criteria were characterized 
using Ingenuity Pathways Analysis (IPA) (Ingenuity Systems, Redwood City, 
CA) functional annotations. Functional annotations were scored using the Fisher 
exact test. Subsequent analysis was performed to better understand the differ¬ 
ences between strains. Textbook ANOVA was performed on the ratio data (i.e., 
ratios determined with respect to the averages of the strain-matched mock- 
infected sample results) at 2 dpi. The following comparisons of mouse groups 
were performed: WT versus IFNAR1“~, WT versus ST ATI -7- , and 
IFNARl -7- versus STAT1 -7- . As described above, a significance criterion of 
P < 0.05 was applied. The gene lists were trimmed by requiring that the change 
of at least one ratio comparing mock results was greater than 1.5-fold. The gene 
lists were then split based on the sample that exhibited the greater magnitude of 
change, yielding 6 gene lists based on both ANOVA and the relative differential 
regulation results for the mouse groups, as follows: genes that are differentially 
regulated more in WT than in IFNARl -7 ', more in IFNARl -/ ~ than in WT, 
more in WT than in STAT1 -7 ”, and so forth. As described above, these gene lists 
were functionally annotated using IPA. This process was repeated for data 
obtained at 5 dpi. The resulting top 5 statistically significant canonical pathways 
for each gene list are shown (see Table 2). Additional analysis of interferon- 
stimulated genes (ISGs) was performed using a panel of genes induced in the 
lungs of WT mice upon intranasal stimulation with 10,000 units of interferon 
alpha for a 24-h time period (14). Gene lists characteristic of T H 1 and T H 2 
immune responses were derived from annotations from the SABiosciences and 
NCBI Entrez Gene databases. Gene lists characteristic of Ml and M2 activated 
macrophages were derived from the work of Mantovani and coworkers (45) and 
NCBI Entrez Gene. 

Masson’s trichrome staining. Lung tissues were fixed in PBS-4% paraformal¬ 
dehyde (pH 7.3) and embedded in paraffin, and 5-|xm-thick sections were pre¬ 
pared by the UNC histopathology core facility. Prior to use, sections were 
deparaffinized in xylene with 2 changes of 5 min each and washed with water. 
Slides were left to incubate overnight in picric acid at room temperature. Slides 
were then washed in running tap water for 10 min to remove the yellow color and 
rinsed with deionized water. The slides were then incubated in Weigert’s iron 
hematoxylin solution (catalog number HT079; Sigma) for 10 min. After incuba¬ 
tion, slides were washed in running water for 10 min and rinsed in deionized 
water. The slides were then incubated in Biebrich scarlet-acid fuschin solution 
(catalog number HT151; Sigma) for 5 min and washed in deionized water. After 
the washing was performed, the slides were incubated in phosphomolybdic- 
phosphotungstic acid (catalog number F9252; Sigma) for 10 min and then in 
aniline blue for 5 min (catalog number 415049; Sigma) before being rinsed in 
deionized water. The slides were then fixed in 1% acetic acid for 3 min followed 
by two rinses of 95% and 100% ethanol and a final xylene rinse. 

Immunohistochemistry. Lung tissues were fixed in PBS-4% paraformalde¬ 
hyde (pH 7.3) and embedded in paraffin, and 5-|xm-thick sections were prepared 
by the University of Maryland histopathology core facility. Sections were depar¬ 
affinized in xylene, with 2 washes of 5 min each. The slides were then placed in 
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FIG. 1. Viral replication in the lungs of wild-type and knockout 
mice. Lung viral titer data from 5 mice infected with rMA15 per strain 
per time point are shown. The limit of detection for this assay is 5 X 
10 2 PFU. 


absolute alcohol and in 95% and 70% alcohol, with 2 washes of 2 min each. After 
2 washes in PBS, the slides were immersed in methanol-0.3% H 2 0 2 for 30 min 
to quench the endogenous peroxidase. The slides were then preincubated with 
10% goat serum for 20 min at room temperature followed by incubation with 
primary antibody for 1 h at room temperature. Rabbit anti-arginase I (Santa 
Cruz Biotechnology, Santa Cruz, CA) was used at a 1:50 dilution, polyclonal 
rabbit anti-YMl (StemCell Technologies, Vancouver, BC, Canada) was used at 
a 1:200 dilution, and monoclonal rat anti-Mac3 (Cedarlane Laboratories, Burl¬ 
ington, NC) was used at a 1:500 dilution. Slides were washed 3 times in PBS, and 
biotin secondary antibody was applied for 30 min at room temperature. The 
sections were washed in PBS 3 times and incubated with ABC reagent (Vector 
Laboratories Inc.). After 2 changes of PBS, DAB (3,3-diaminobenzidine) was 
used to develop the sections, which were then counterstained with Mayer’s 
hematoxylin for 2 min. Slides were dehydrated in 2 changes of 95% and absolute 
alcohol and then mounted with coverslips. 

Quantitative RT-PCR. Relative quantitative reverse transcription-PCR (RT- 
PCR) was performed to validate a number of gene expression changes as de¬ 
tected with microarrays. A QuantiTect reverse transcription kit (Qiagen Inc., 
Valencia, CA) was used to generate cDNA, and quantitative RT-PCR was 
performed using an Applied Biosystems 7500 PCR system. The following 
TaqMan assays for mouse genes were obtained from Applied Biosystems: for 
Igfl, assay Mm00439560_ml; for Cclll, assay Mm00441238_ml; for Argl, assay 
Mm00475988_ml; for 1117a, assay Mm00439619_ml; and for Col3al, assay Mm 
01254476_ml. The TaqMan assay for SARS nucleoprotein (NP) was described 
previously (39). Loading variations were normalized using the invariant calibra¬ 
tor gene Mfap-la according to the 2 _AACT threshold cycle method (41). Differ¬ 
ences in gene expression are represented as log ratios relative to cDNA samples 
from strain-matched mock-infected animals. 

FACS analysis of alternatively activated macrophages (AAMs). Ten-week-old 
129SvEv and STATU - mice (Taconic Laboratories) were infected with 10 5 
PFU of rMA15 as described previously. Mice were euthanized at day 8 postin¬ 
fection and their lungs dissected for fluorescence-activated cell sorter (FACS) 
analysis. Lungs of the mice were dissociated in collagenase and DNase, after 
which the cells were incubated with PerCP-Cy5.5-labeled anti-CDllb (eBio- 
science, San Diego, CA) and phycoerythrin-labeled anti-CD 14 (eBioscience, San 
Diego, CA). Cells were permeabilized and stained with anti-Mac3 (Cedarlane 
Laboratories Ltd., Burlington, NC) and then with allophycocyanin-labeled goat 
anti-rat IgG (Biolegend, San Diego, CA). Samples were processed on a Becton 
Dickinson LSR II FACS machine and analyzed using FlowJo 7.6. 


RESULTS 

To better understand the contributions of specific signaling 
pathways and gene expression changes to SARS-CoV-induced 
lung disease, we examined the global transcriptional changes 
that occur during rMA15 infection of WT mice or mice con¬ 
taining genetic defects associated with innate immune signal¬ 
ing. WT, IFNAR1 and STAT1 “ mice were infected with 
1 X 10 5 PFU of rMA15. Five animals per group were sacrificed 


at 2, 5, and 9 dpi. Weight loss, tissue pathology, and viral titer 
data are shown in Fig. 1 and 2 and Table 1; a detailed descrip¬ 
tion of the disease course was previously reported (24). Briefly, 
WT and IFNARl 7- mice experienced similar disease courses 
based on viral titer, weight loss, and lung pathology results. WT 
and IFNAR1 “ mice lost 15% of their starting weight by 5 dpi 
but regained this weight by 9 dpi. Peak viral titers in the lung 
for both mouse strains were observed at 2 dpi and diminished 
over the course of the infection. At 9 dpi, virus was not de¬ 
tectable by plaque assay or by in situ hybridization. 

STATU 7- mice exhibited weight loss and pathology similar 
to the results seen with the WT and IFNARl -7- mice at 2 dpi 
but experienced increased weight loss and inflammatory cell 
infiltration into the lungs at 5 and 9 dpi. Notably, extensive 
collagen deposition was observed in the lungs of STAT1 ~ 
mice by Frieman et al. (24) and in this work as judged by 
Masson’s trichrome staining (Fig. 3), suggesting early-stage 
pulmonary fibrosis. STAT1 -7- mice also exhibited higher viral 
titers in the lung throughout the time course, including 5 X 10 6 
PFU per gram of tissue at 9 dpi. Gene expression profiling 
analyses of RNA samples from 3 of the 5 mice per group were 
used to identify mechanisms of protection against rMA15 ex¬ 
hibited by wild-type and IFNARl -7- mice and to identify the 
pathological mechanisms that lead to fibrosis-like disease in 
the STATU 7- mice, as described in the following sections. 

Immune and mitotic gene expression reveals common anti¬ 
viral response and damage repair mechanisms. We began our 
gene expression studies by performing ANOVA of gene ex¬ 
pression data obtained from the lungs of WT, IFNARl -7- , 
and STAT1 -7- mice. ANOVA identifies genes for which the 
variation in expression levels between biological replicates is 
significantly smaller than the variation between different treat¬ 
ment groups (based on either time or mouse strain). We initially 
focused on identifying functional gene classes that were associ¬ 
ated with the progression of rMA15 infection and pathology. To 
accomplish this, we performed three comparisons for each strain: 
(i) expression data from mock-infected animals were compared 
with expression data from infected animals sacrificed at 2 dpi; (ii) 
expression data from infected animals sacrificed at 2 dpi were 
compared with expression data from infected animals sacrificed at 
5 dpi; and (iii) expression data from infected animals sacrificed at 
5 dpi were compared with expression data from infected animals 
sacrificed at 9 dpi. 

Using IPA to evaluate the gene sets determined using our 
ANOVA results, we identified responses common to rMA15- 
infected WT, IFNARl -7- , and STAT1 -7- mice relative to 
mock-infected animals. The commonly identified response at 2 
dpi consisted of a set of genes involved in immunity and the 
inflammatory response to infection that were annotated by the 
IPA “immunological disease” functional category. WT mice 
differentially regulated 183 genes (P = 1.85 X 10 -36 ), 
IFNARl 7- mice differentially regulated 230 genes (P = 
4.36 X 10“ n ), and STAT1 -7- mice differentially regulated 202 
genes (P = 8.30 X 10 45 ) in this category. A union of these 
genes is shown as a heat map in Fig. 4 (top panel). Increased 
expression of a subset of these genes at 2 dpi, including genes 
encoding cytokines (Ilia, Illb, 116, 1110, 1112b, and Tnf), che- 
mokines (Ccl2, Ccl3, Ccl4, Ccl7, Ccl20, Cxcl2, Cxcl3, and 
CxcllO), and immune sensor molecules (Tlrl, Tlr2, and Tlr3), 
suggested a common initial immune and inflammatory re- 
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FIG. 2. Sustained inflammation in the lungs of STAT1 7 mice infected with rMA15. Representative tissue samples from WT, IFNAR 7 A and 
STAT1 7 mice infected with rMA15 or mock-infected (2, 5, and 9 dpi) were stained with hematoxylin and an eosin counterstain. Images at lower 
magnification are available in Fig. SI in the supplemental material. Magnification, X40. 


sponse in all three strains of mice (Fig. 4, top panel [yellow 
box]) (individual genes are listed in Table SI in the supple¬ 
mental material). 

By 5 dpi, the mRNA abundance of these same immune- 
related genes was reduced in WT and IFNAR1 /_ mice, which 
correlated with the reduced pathology seen with these animals. 
In contrast, expression of this gene set in STAT1 “ mice was 
increased at 5 dpi, reflecting an observed increase in inflam¬ 


mation. At 9 dpi, expression of immune-related genes was 
further reduced in WT and IFNAR1 “ mice whereas expres¬ 
sion of these genes was still significantly increased in 
STAT1 -- mice. This may reflect the extensive immune cell 
infiltration observed in the lungs of these animals (24). 

Our ANOVA and IPA results also suggested the presence of 
extensive cell proliferation by 5 dpi. In particular, the IPA “cell 
stage” functional category was strongly associated with this 


TABLE 1. Summary of clinical and virological measures 


Characteristic 


Value for indicated mouse strain 13 


129 

IFNAR1” 7 ” 

STAT1 -7- 

Time of maximum weight loss (dpi) 

4 

4 

9 

Maximum wt loss (% of starting wt) 

14-16 

14-16 

>20 

Peak viral titer (PFU/g) 

2.6-3.8 x 10 7 

3-7.2 x 10 7 

7.6 X 10 7 -3.4 x 10 8 

Viral RNA level at 2 dpi 

6.4-16 x 10 4 

1.4-2.8 x 10 5 

8.9-23 x 10 4 

Viral RNA level at 5 dpi 

1.2-20 x 10 2 

3.7-25 x 10 2 

3.1-4.4 x 10 4 

Viral RNA level at 9 dpi 

1.0-7.5 

5.2-84 

5.5-9.5 x 10 3 

Alveolar protein deposition at 9 dpi 

No 

No 

Yes 

Endpoint 

Recovery 

Recovery 

Death 


a The ranges of viral RNA levels are shown on an arbitrary scale as determined by quantitative RT-PCR results determined using probes for the SARS nucleoprotein 
gene (39). Three animals were used per set of conditions and time point. 
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STAT1 -/- + PBS STAT1 -/- + rMAI 5 



FIG. 3. Abnormal protein deposition in the lungs of STAT1 -/ 
mice infected with rMA15. Representative tissue samples from 
STAT1 -/ mice infected with rMA15 or subjected to mock infection (9 
dpi) were stained with Masson’s Trichrome. The presence of blue- 
purple staining indicates protein deposition in the alveolar spaces. 
Images at lower magnification are available in Fig. S2 in the supple¬ 
mental material. 


time point. WT mice differentially regulated 171 genes (P = 
1.05 X 10 _3 °), IFNAR1 “ mice differentially regulated 185 
genes (P = 1.82 X 10“ 19 ), and STAT1 -- mice differentially 
regulated 160 genes (P = 4.71 X 10“ 14 ) in this category. 
Expression of a union of these genes is shown in Fig. 4 (lower 
panel). Of particular interest are members of a subset of these 
genes that exhibit low expression at 2 dpi and high expression 
at 5 dpi, with the results being independent of the mouse strain 
(Fig. 4, cyan box). These genes include DNA damage repair- 
associated genes (Brcal, Rad51, and Rad541), cyclins (Ccna2, 
Ccnbl, Ccnel, Ccne2, and Ccnf), cell division cycle genes 
(Cdc6, Cdc20, Cdc25c, Cdca5, and Cdca8), and genes associ¬ 
ated with DNA replication initiation (McmlO and Pole). A full 
list of these genes is included in Table S2 in the supplemental 
material. In particular, IFNAR1 /_ mice exhibited the stron¬ 
gest expression of these genes, suggesting that more cell divi¬ 
sion occurs in these animals at 5 dpi. By 9 dpi, the expression 
of these genes was reduced in WT and IFNAR1 /_ mice, 
suggesting a decrease in cell proliferation as the virus was 
cleared. In contrast, STAT1 /_ mice exhibited a further in¬ 
crease in gene expression associated with the cell cycle at 9 dpi. 
This may have been due to proliferation of fibroblasts or cell 
division of infiltrating immune cells. 

Differential host response patterns suggest mechanisms of 
pathogenesis and survival. In order to reveal the differences in 
the host response that underlie the inability of STAT1 /_ mice 
to clear the virus and the subsequent development of pulmo¬ 
nary fibrosis-like lesions, we performed ANOVA to identify 
genes that were differentially expressed based on the mouse 
strain. These genes were then functionally annotated using 
pathway information from IP A. Table 2 contains the pathways 
most significantly enriched in the ANOVA gene lists from 2 
dpi. Genes associated with the interferon signaling pathway, 
including Ifi35, Ifitl, Ifit3, Irfl, Irf9, Jak2, Medl4, Oasl, 
Psmb8, Socs3, Statl, Stat2, and Tapi, were robustly upregu- 




Log(Ratio) 

- 1.0 


FIG. 4. Kinetics of the transcriptional response to infection corre¬ 
late with pathological findings. Expression profiles of 261 immune- 
related genes (top) and 204 cell cycle-related genes (bottom) selected 
by ANOVA are shown. Each column represents averages of data from 
the results of 3 replicate animal experiments performed using WT, 
IFNAR1 _/_ , or STAT1 7 mice at 2, 5, or 9 dpi compared to averages 
of data from experiments performed using mock-infected mice of the 
same strain. Genes shown in red were upregulated, genes shown in 
green were downregulated, and genes shown in black indicate no 
change in expression relative to mock-infected mice. Genes in the top 
panel were identified as part of the IPA ‘‘immunological disease” 
functional category in at least one ANOVA comparison between 
mock-infected animals and animals sacrificed at 2 dpi. A yellow box 
highlights 81 genes similarly upregulated by the 3 strains at 2 dpi. The 
identities of these genes are listed in Table SI in the supplemental 
material. Genes in the lower panel were identified as part of the IPA 
“cell stage” functional category in at least one ANOVA comparison 
between mice sacrificed at 2 and at 5 dpi. A cyan box highlights 80 
genes commonly upregulated by all strains of mice by 5 dpi. The 
identities of these genes are listed in Table S2 in the supplemental 
material. 


lated in WT or IFNAR1 /_ mice compared with STAT1 
mice (Table 2). 

Since most of the genes noted above are induced in response 
to interferon stimulation, we compared our ANOVA gene list 
with a more comprehensive list of genes that have been shown 
to be upregulated in mice treated intranasally with interferon 
alpha (see Materials and Methods and reference 14). We iden- 
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TABLE 2. Top differentially expressed pathways at 2 dpi 


Gene regulation comparison for 
indicated strains" 

Canonical pathway 6 

No. of 
genes 

Score" 

129 more regulated than IFNAR1 7 

Acute-phase response signaling 

15 

4.3 


LXR/RXR activation 

9 

4.1 


FXR/RXR activation 

10 

4.0 


TREM1 signaling 

8 

4.0 


Hepatic cholestasis 

11 

3.1 

129 more regulated than STAT1 7 

Pattern recognition receptors in recognition of bacteria and viruses 

21 

7.6 


Interferon signaling 

11 

5.9 


Allograft rejection signaling 

12 

5.3 


Antigen presentation pathway 

11 

5.3 


Graft-vs-host disease signaling 

12 

5.2 

IFNAR1 7 more regulated than 129 

Valine, leucine, and isoleucine degradation 

25 

7.4 


Propanoate metabolism 

21 

5.9 


LPS/IL-1-mediated inhibition of RXR function 

43 

4.9 


(3-alanine metabolism 

18 

4.8 


Pyruvate metabolism 

22 

4.8 

IFNAR1 7 more regulated than STAT1 -7 

Valine, leucine, and isoleucine degradation 

28 

8.9 


Propanoate metabolism 

23 

6.1 


Pyruvate metabolism 

23 

5.5 


(3-alanine metabolism 

18 

4.4 


Interferon signaling 

12 

4.1 

STAT1 -/ more regulated than 129 

Arginine and proline metabolism 

9 

4.3 


Urea cycle and metabolism of amino groups 

6 

4.2 


Bile acid biosynthesis 

7 

4.0 


IL-10 signaling 

7 

3.2 


Complement system 

5 

2.9 

STAT1 7 more regulated than IFNAR1 7 

LXR/RXR activation 

9 

6.2 


Acute-phase response signaling 

12 

5.2 


IL-10 signaling 

7 

4.6 


Cytokine mediation of communication between immune cells 

6 

3.9 


Hepatic cholestasis 

9 

3.8 


a Lists of differentially expressed genes were determined as described in Materials and Methods and were investigated using IPA. The 5 top-scoring canonical 
pathways are shown. 

b LXR, liver X receptor; RXR, retinoid X receptor; FXR, farsenoid X receptor; LPS, lipopolysaccharide. 

c Enrichment scores are calculated as —logfP value). The P value is derived from Fisher’s exact test as described in Materials and Methods. 


tified 268 interferon-stimulated genes that were more strongly 
induced or repressed in WT or IFNAR1 /_ mice than in 
STAT1 -7- mice (Fig. 5, orange and yellow bars). Genes an¬ 
notated with an antiviral function were most significantly en¬ 
riched in this gene set (Fig. 5, yellow bar), suggesting that 
signaling independent of the type I interferon receptor, but 
dependent on STAT1, activates interferon-stimulated gene 
(ISG) expression during rMA15 infection, including expression 
of genes required for clearance of rMA15. Such signaling could 
be the result of interferon lambda production. 1128b (inter¬ 
feron lambda 3) is modestly induced at 2 dpi, averaging 4.5- 
fold induction in 129 mice, 7.4-fold induction in IFNAR1 7 
mice, and 10.2-fold induction in STAT1 7_ mice but not at 
later time points (data not shown). This level of induction is 
modest compared with the more than 10-fold induction seen 
with a number of type I interferons at 2 dpi (data not shown). 

In contrast, only 48 ISGs were identified as being more 
strongly induced in STAT1 7 mice than in WT or 
IFNAR1 “ mice (Fig. 5, blue and cyan bars). Many of these 
genes are involved in leukocyte trafficking (Fig. 5, cyan bar). In 
particular, S110a9, Csf3r, Mmp8, Cxcll, Saal, CxcllO, Ccl7, 
and 116 gene expression is associated with neutrophil infiltra¬ 


tion (4, 10, 22, 48, 55, 74), Ccl8 and Cclll are associated with 
eosinophilia (10), and Ccl4 is associated with macrophage ac¬ 
tivation (20). Notably, levels of neutrophils, eosinophils, and 
macrophages were elevated in the lungs of STAT1 7_ mice 
infected with rMA15 (24). These results suggest that leukocyte 
infiltration of the lung is the main antiviral response available 
to STAT1 “ mice in the absence of effective ISG induction. 
Although these genes were identified as stimulated by inter¬ 
feron, many of these genes are also activated by other proin- 
flammatory signaling pathways, most notably AP-1 and NFkB 
(18, 54, 72). 

STAT1 -7- mice exhibit induction of genes associated with 
leukocyte recruitment, T H 2 bias, and fibrosis. In order to un¬ 
derstand the regulatory patterns associated with the increased 
pathogenesis of rMA15 in STAT1 -7- mice, we examined those 
genes selected by ANOVA with the STAT1 7_ and WT or 
IFNAR1 “ datasets (described above) that exhibited greater 
differential regulation in STAT1 7 “ mice. At 2 dpi, this anal¬ 
ysis revealed a significant enrichment for genes associated with 
leukocyte trafficking (Table 2). In addition, eosinophilic che- 
mokine Ccl24 and T H 2 signaling transcripts (Csf3,115, and 119) 
were identified. The involvement of T H 2-associated genes was 
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FIG. 5. Deviant interferon response in STAT 1 7 mice. Expression 
profiles are shown for interferon-stimulated genes (identified by intrana¬ 
sal administration of 10,000 units of interferon alpha for 24 h) (see details 
in Materials and Methods and reference 14) that were significantly dif¬ 
ferent between STAT1 7 mice and WT or IFNAR1 7 mice as judged 
by ANOVA. Each column represents gene expression data from an av¬ 
erage of 3 animal replicate experiments compared to mock-infected ani¬ 
mal results. Genes shown in red were upregulated, genes shown in green 
were downregulated, and genes shown in black indicate no change in 
expression relative to mock-infected animal results. Genes are marked 
according to whether greater dilferential regulation was observed in 
STAT1 7 mice (blue and cyan) or in WT and IFNAR1 7 mice (orange 
and yellow). Genes marked with a cyan bar were functionally annotated as 
involved with “cell movement of leukocytes” by IPA (14 genes; P = 
8.61 X 10 1 '). Genes marked with a yellow bar were functionally anno¬ 
tated as “antiviral response” genes by IPA (18 genes; P = 8.58 X 10 12 ). 
Genes in these functionally annotated categories are listed by Entrez 
Gene name. Genes comprising this heat map are listed in Table S3 in the 
supplemental material. 


unexpected in the context of an acute viral infection of naive 
animals. T H 2 responses are associated with humoral immunity 
to bacteria and viruses, wound healing, and helminth infection, 
whereas clinical and experimental SARS infections are gener¬ 
ally characterized by a T H 1 immune response (40, 63). In 
addition, STAT1 /_ mice exhibited increased gene expression 
associated with the arginine-proline metabolism pathway at 2 
dpi (Table 2), including inducible arginase (Argl), a marker 
for profibrotic “alternatively activated macrophages,” as de¬ 
scribed in detail below. 

Of particular significance, numerous genes associated with 
the development of fibrosis were differentially induced by 


WT IFNAR1-/- STAT 1 -/- 
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FIG. 6. Increased fibrotic gene expression in STAT1 7 mice. Ex¬ 
pression profiles are shown for 48 genes associated with the IPA 
“hepatic fibrosis” canonical pathway. The heat map is colored accord¬ 
ing to the ratio of gene expression at 5 dpi to that seen with mock- 
infected mice. Red indicates upregulation, green indicates downregu- 
lation, and black indicates no change in expression upon infection. 


STAT1 “ mice compared with WT and IFNAR1 “ mice at 
5 dpi (Fig. 6), including mediators of inflammation (Ccl2, 
Cxcl3, Ifng, Ilia, Illb, 116, Nfkb2, and Tnf) and tissue remod¬ 
eling (Igfl, Igfbp4, Igfbp3, Timpl, Timp2, Mmp9, Mmpl3, 
Fnl, and Col3al). These results suggested that we were able to 
detect the activation of profibrotic gene expression several days 
prior to observing a fibrotic phenotype in the infected lung; we 
previously observed a similar phenomenon in the context of he¬ 
patic fibrosis caused by hepatitis C virus (68). RT-PCR validation 
performed using Col3al and Igfl confirmed this pattern of gene 
expression (data not shown). Fibrosis genes included extracellular 
matrix proteins (Col3al and Fnl), metalloprotease inhibitor 
Timpl, and insulin-like growth factor (Igfl), which have previ¬ 
ously been implicated as hallmarks of pulmonary fibrosis progres¬ 
sion in the bleomycin model (11, 43, 70). 

Knockout mice exhibit a bias in T-helper and macrophage 
differentiation. Because several genes identified in the time- 
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FIG. 7. Increased T H 2- and M2-associated gene expression in 
STAT1 7 mice. (A and B) Expression profiles are shown for genes 
associated with T H 1 and T H 2 responses (A) and for Ml classical 
macrophage activation and M2 alternative macrophage activation (45) 
(B). Individual genes are identified by the mouse Entrez Gene name. 
T h 2 responses and M2 responses were upregulated only in STAT1 7 
mice and are highlighted by yellow and cyan boxes, respectively. 


matched ANOVA were associated with a T H 2 immune re¬ 
sponse, we investigated gene regulation associated with the 
T H 1/T H 2 balance in greater detail. Expression levels for T H 1- 
associated genes (Fig. 7A) were similar in all mouse strains at 
2 dpi. However, the STAT1 7_ mice exhibited stronger induc¬ 
tion of T h 1 genes at 5 and 9 dpi compared with WT and 
IFNAR1 7_ mice. This may reflect the ongoing viral replica¬ 
tion and activation of the T H 1 master regulator T-bet in the 
lungs of STAT1 “ mice. In contrast, whereas WT and 
IFNAR1 7_ mice exhibited little regulation of T H 2-associated 
gene expression, STAT1 “ mice exhibited significant induc¬ 
tion of several T H 2 genes, including 115,119, Ccl7, Cclll, Ccl24, 
and Csf3 (Fig. 7A, yellow box). RT-PCR performed on Cclll 
confirmed this pattern of gene regulation (data not shown). 
These data suggest that STAT1 7_ animals experience a T H 2 
bias compared with WT mice. 

Because STAT1 7_ mice exhibited dysregulation of the T H 1 
and T H 2 responses, we examined the possibility that the re¬ 
cently described T H 17 lineage could also be perturbed in these 
animals. These cells contribute to pathogenesis in asthma and 
chronic inflammation (49, 73, 81) and could be relevant to the 
ongoing pathogenesis observed in STAT1 -7- mice. Although 


induction of T H 17 cytokine gene 1117a was negligible in WT 
and IFNAR1 7_ mice, expression of this gene was strongly 
induced at 9 dpi in STAT1 7 “ mice (averaging 86-fold by 
microarray and 1,600-fold by RT-PCR) (data not shown). 
Likewise, the expression of 1123a, a cytokine responsible for 
survival of T H 17 cells, was induced 4-fold (based on microarray 
data) in STAT1 7_ animals at 9 dpi. These results suggest that 
activation of the T H 17 pathway occurs concurrently with the 
development of fibrosis in STAT1 7 “ mice. 

Effect of the T H 2 bias on macrophage phenotype. Previous 
studies have suggested that macrophage phenotypes play a role 
in T H 2-mediated fibrosis (2, 25, 51). In response to viral or 
bacterial pathogens, macrophages become activated by gamma 
interferon along a classical pathway (Ml macrophages). These 
cells are characterized by production of proinflammatory cy¬ 
tokines as well as nitric oxide synthase (45). Gene expression 
for individual Ml-associated genes was strong at 2 dpi in all 
mouse strains studied (Fig. 7B). This gene expression de¬ 
creased in the WT and IFNAR1 “ mice at 5 and 9 dpi as 
these strains cleared the virus. In contrast, Ml gene expression 
remained elevated in STAT1 7_ mice, reflecting sustained in¬ 
fection in these animals. 

Under certain circumstances, such as response to helminth 
infections, macrophages (specifically, M2 macrophages) can 
become activated by 114 along an alternative pathway. These 
AAMs are characterized by production of arginase I, 1110, 
eosinophilic chemokines, and insulin-like growth factor-1. In¬ 
creased expression of individual M2-associated genes, includ¬ 
ing arginase I, Ccl24, Cdl4, and I18rb by 2 dpi and Cell, Igfl, 
and Ccr2 by 5 dpi, revealed strong activation of the M2 re¬ 
sponse only in STAT1 7_ animals (Fig. 7B, cyan box). This 
pattern of gene expression was also confirmed using RT-PCR 
for Argl and Igfl (data not shown). 

Lungs of STAT1 -7- mice exhibit macrophage infiltration 
and activation along the alternative pathway. To confirm the 
production of arginase I by M2 macrophages in STAT1 -7- 
mice, we performed immunohistochemistry to detect several 
protein markers of alternative macrophage activation: arginase 
I, Yml, and Mac3. As shown in Fig. 8, minimal staining was 
detected in the lungs of mock-infected animals. By 9 dpi, both the 
number of cells staining for these proteins and the density of the 
staining had increased in STAT1 mice. Increased staining was 
also observed at 2 and 5 dpi (data not shown). Of particular 
interest, the cells exhibiting high levels of staining exhibited mor¬ 
phology characteristics consistent with macrophages. These re¬ 
sults extend our gene expression data, suggesting that, in 
STAT1 7_ mice, arginase I-expressing M2 macrophages infiltrate 
the lung in large numbers during rMA15 infection. 

To confirm that lung macrophages were being activated 
along the alternative pathway, we performed flow cytometry 
using antibodies to CD lib and CD 14 to identify macrophages 
and Mac3 to identify AAMs (A. Q. Ford, P. Dasgupta, E. A. 
Smith, N. Noben-Trauth, and A. D. Keegan, submitted for 
publication). Numbers of AAMs remained low in strain 129 
mice and uninfected STAT1 -7- mice. However, STAT1 -7- 
mice infected with rMA15 exhibited a dramatic increase in the 
population of AAMs (Fig. 9A) as well as in the number of 
AAMs as a percentage of total macrophages (Fig. 9B). Based 
on these data, we conclude that AAMs are disproportionately 
activated in STAT1 -7- mice infected with rMA15. 
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FIG. 8. Detection of AAM markers in lung tissue of infected 
STAT1 -/ mice by the use of immunohistochemistry. Lung sections 
were stained for the AAM markers arginase I, Yml, and Mac3 with a 
hematoxylin counterstain. At 9 dpi, prominent staining is observed in 
STAT1 7 mice infected with rMA15 compared with uninfected ani¬ 
mals and WT animals infected with rMA15. Magnification, X40. 


DISCUSSION 

Previous studies have demonstrated that SARS-CoV infec¬ 
tion of mice mirrors the acute exudative phase of ARDS (63, 
64). More recently, STAT1 7_ mice were shown to be highly 
susceptible to rMA15 infection, whereas disruption of inter¬ 
feron signaling at the level of the type I, II, or III interferon 
receptors did not significantly alter the disease course from 
that observed in WT mice (24). Our current gene expression 
analyses further demonstrate that although there were differ¬ 
ences in the transcriptional profiles of WT and IFNAR1 “ 
mice, signaling through the type I interferon receptor was not 
necessary for a protective antiviral transcriptional response 
during rMA15 infection. Other studies performed by our 
group have revealed that IFNAR1 7_ mice, and mouse em¬ 
bryonic fibroblasts derived from these animals, also induce ISG 
expression in response to infection with influenza virus (14, 
27), suggesting that there exists a common mechanism of com¬ 
pensation for the loss of type I interferon receptor signaling. 
Our results also reveal molecular mechanisms by which 


A 



PBS MAI 5 PBS MAI 5 
129 STAT1 -/- 



PBS MAI 5 PBS MAI 5 
129 STAT 17- 

FIG. 9. Enrichment of alternatively activated macrophages in 
STAT1 7 mice infected with SARS-CoV. Ten-week-old strain 129 
WT and STAT1 7 mice were infected with MA15, and their induction 
of AAMs was analyzed at 8 days postinfection. A total of 50,000 lung 
cells were sorted for the total macrophage population according to the 
presence of CD lib and CD 14 double positives. Of these cells, AAMs 
were identified by the presence of Mac3 in the cells by intracellular 
staining. (A) The numbers of CDllb + , CD14 + , and MAC3 + cells are 
graphed. (B) The numbers of AAMs are graphed as percentages of the 
total number of macrophages (CDllb+/CD14 + cells). 


STAT1 7_ mice progress to fibrosis-like disease, mirroring the 
exudative and fibrotic phases of ARDS. 

Compensatory gene expression in IFNAR1 - - mice. Al¬ 
though WT and IFNAR1 7_ mice exhibited similar weight loss 
and pathology characteristics, absence of the type I interferon 
receptor was associated with changes in host gene expression 
with respect to the immune response and healing. For exam¬ 
ple, at 2 dpi, a large number of genes associated with metab¬ 
olism exhibited decreased expression in IFNAR1 7_ mice (Ta¬ 
ble 2). Similar metabolic changes have been observed in 
response to a variety of viral infections, suggesting that this is 
an additional protective response to viral infection (3, 8, 13). 
Although genes associated with mitosis were induced in WT 
and IFNAR1 7_ mice at 5 dpi, IFNAR1 7 “ mice exhibited 
greater induction of these genes, particularly those associated 
with DNA damage repair. In addition, we observed a slight 
attenuation of ISG expression in IFNAR1 7_ mice compared 
with WT mice (Fig. 5). This attenuation may be responsible for 
the larger viral load observed in the lungs of IFNAR1 7 “ mice 
at 5 dpi compared with that observed in WT mice at the same 
time point (24). This compensatory ISG activation was poten¬ 
tially the result of interferon lambda signaling, but only modest 
induction of 1128b (interferon lambda 3) was observed in the 
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FIG. 10. Dysregulation of T-helper and macrophage differentiation 
in STAT1 7 mice. A functional network that depicts the effects of 
T-helper differentiation on macrophage phenotype and disease reso¬ 
lution is shown. The regulation of genes associated with T H 1, T H 2, Ml, 
and M2 is represented with red arrows to indicate upregulation or with 
a gray bar to indicate no significant change from mock-infected animal 
results. A heat map is shown for select macrophage-associated genes 
that contribute to fibrosis and eosinophilia. Genes are colored accord¬ 
ing to differential regulation at 9 dpi. 


lungs of all strains at 2 dpi. However, Frieman et al. showed 
that mice with defects in both type I and type III interferon 
signaling cleared rMA15 infection, suggesting the presence of 
additional compensatory pathways (24). We also cannot rule 
out the possibility that induction of some ISGs was the result of 
secondary signaling through other pathways, including that 
involving interferon gamma. However, Frieman et al. showed 
that interferon gamma signaling was also dispensable for clear¬ 
ing rMA15 (24). Through redundant pathways, some defects in 
innate immunity can be overcome, including removal of type I 
interferon signaling. By comparison, there was a dramatic re¬ 
duction of ISG expression in rMA15-infected STAT1 -7- mice; 
this reduction correlated with increased pathology and death 
for these animals. Since many immune response pathways, 
including the interferon response, signal through STAT1, the 
STAT1 7 “ mice appear to be unable to mount an effective 
antiviral response to rMA15 infection. 

Macrophage dysregulation in STAT1 - - mice. Our gene 
expression data suggest that infection of STAT1 7_ mice with 
rMA15 dysregulates T H 1/T H 2 and macrophage host responses, 
resulting in a profibrotic phenotype. A model consistent with 
our data is shown in Fig. 10. In WT and IFNAR1 “ mice (but 
not STAT1 7 “ mice), T cells that differentiate with a T H 1 
program repress T H 2 differentiation via gamma interferon sig¬ 
naling through STAT1 (82). Gamma interferon also stimulates 
the differentiation of macrophages along a classical activation 
pathway (Ml) characterized by production of interleukins 1 
and 6, tumor necrosis factor, and Cxcr3 ligands (Cxcl9, CxcllO, 
and Cxclll) (45); upregulation of these genes was observed in 
WT and IFNAR1“ 7 ~ mice at 2 dpi. 

In contrast, we found that STAT1mice exhibited a T H 2- 
biased immune response to infection with rMA15. T H 2 cells 


produce 114 (production of which was transcriptionally in¬ 
creased in STAT1 7_ mice infected with rMA15), which trig¬ 
gers the differentiation of macrophages along an alternative 
activation pathway (M2). Alternatively activated macrophages 
were observed at 2 dpi in STAT1 7_ mice, and their numbers 
increased at 5 and 9 dpi. The T H 2/M2 response is associated 
with a variety of biological functions, including wound healing, 
chronic inflammation, and fibrotic immune responses to hel¬ 
minths (15, 66, 78). Notably, similarities with respect to 
T H 2/M2 bias and pulmonary fibrosis have been observed in 
gamma interferon receptor knockout (IFNGR 7_ ) mice in¬ 
fected with murine gammaherpesvirus 68 (51). In contrast, we 
did not observe fibrosis in IFNGR 7_ mice infected with 
rMA15 (24). This could have been due either to inhibition of 
uncontrolled cell proliferation by STAT1 (76) or to effective 
clearance of the virus in these mice. 

Thus, a signaling pathway that promotes healing under fa¬ 
vorable conditions may instead promote fibrosis in the context 
of continuing immune dysregulation, uncontrolled fibroblast 
proliferation, and/or persistent viral replication (Fig. 10). One 
hallmark of M2 macrophages is the production of arginase I 
(31), which regulates the first committed step in the conversion 
of arginine to proline (42). Proline is an important building 
block for collagen production, so increases in proline concen¬ 
trations are considered profibrotic. It is interesting that, al¬ 
though arginase I was strongly induced transcriptionally at all 
time points in STAT1 7_ mice, we observed a lag in the pro¬ 
duction of arginase I protein in the lung (Fig. 8). Likewise, 
genes in a panel associated with fibrosis progression (Fig. 6) 
were differentially induced in STAT1 7_ mice prior to the 
appearance of protein deposition in the lungs. Similar results 
have also been found in the context of hepatic cell apoptosis 
and fibrosis (68). These results suggest the potential for diag¬ 
nostic tests to detect gene expression changes predictive of the 
initiation of fibrosis. 

STAT1 7_ mice also exhibited increased expression of ad¬ 
ditional M2 macrophage-associated genes that contribute to 
fibrosis and eosinophilia. For example, we observed increased 
expression of Igfl, an important regulator of cell proliferation 
(65) that promotes fibroblast survival and increased extracel¬ 
lular matrix protein deposition (79). It is negatively regulated 
by several IGF-binding proteins (Igfbpl through Igfbp6), and 
expression of genes encoding these proteins was decreased in 
STAT1 “ mice in response to rMA15 infection (data not 
shown). In addition, we observed increased expression of 
Mmp9 and Mmpl2, which are expressed by M2 macrophages 
to promote wound healing (56) and tissue remodeling (33), 
respectively. 

Contribution of eosinophils and T H 17 cells. T H 2-mediated 
inflammation causes eosinophil infiltration, which is generally 
not observed during acute viral infection. Mmpl2 is an impor¬ 
tant contributor to airway eosinophilia in mice (60), and eo¬ 
sinophilia is associated with pulmonary fibrosis in humans (21) 
and with ARDS (35) and is also observed in some cases of 
persistent inflammation. Expression of genes encoding eosin¬ 
ophil attractants, including Cclll, Ccl24, and Mmpl2, was 
observed in STAT1 7_ mice as early as 2 dpi, and significant 
eosinophilia was observed at 8 dpi (24). 

T h 17 cells are a recently described T-helper population dis¬ 
tinct from T h 1 and T H 2 cells (57) and are associated with 
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asthma, chronic obstructive pulmonary disease, and autoim¬ 
mune inflammation (30, 36, 53, 73, 81). As with eosinophilia, 
the presence of T H 17 cells, which mediate inflammation 
through the production of 1117 and 1123, is not characteristic of 
most acute viral infections. However, gene expression charac¬ 
teristic of T h 17 cells was observed in STATU 7- mice at 9 dpi. 
As with T h 2, T h 17 differentiation is repressed by gamma in¬ 
terferon signaling through STAT1 (57), suggesting that re¬ 
moval of this block in STATU 7- mice enables the expansion 
of T h 17 cells. The occurrence of eosinophils and T H 17 cell 
markers in the lungs of STAT1 -7- mice at 8 to 9 dpi suggests 
that these cells may play a role in the later progression rather 
than the initiation of fibrosis, although this has not been di¬ 
rectly examined. 

Comparison of rMA15-infected STAT1 -7- mice with other 
SARS-CoV models. This is the first report of a study examining 
the host transcriptional responses elicited by the mouse- 
adapted rMA15 variant of SARS-CoV in young mice. The 
disease course caused by rMA15 infection in young strain 129 
mice was similar to that observed in aged BALB/c mice in¬ 
fected with SARS-CoV Urbani (63). Of particular interest, a 
comparison of T H 2 and M2 gene expression characteristics in 
aged and young BALB/c mice (63) reveals greater induction of 
I113ral, Ccl5, Ccl7, Fcerlg, and Igfl in aged mice, suggesting 
that T H 2 and M2 induction upon viral infection is a common 
feature of aged BALB/c and STAT1 7- mice. In contrast, 
induction of T H 17-associated genes was not observed in aged 
BALB/c mice infected with SARS-CoV. In addition, previously 
described SARS-CoV mouse models that used stain rMA15 or 
Urbani (62-64, 67) exhibited either mortality or recovery from 
infection within 6 dpi. In contrast, STAT1 7- mice did not 
succumb to the infection until approximately 10 dpi and ex¬ 
hibited protein deposition in the lungs consistent with fibrosis. 
This suggests that rMA15-infected STAT1 -7- mice represent a 
model of both the exudative and the organizing (fibrotic) 
phases of ARDS pathophysiology (61). 

Previous transcriptome analysis of blood samples from pa¬ 
tients suffering from SARS-CoV infection suggested that sus¬ 
tained innate immune activation and a failure to switch to 
humoral immunity are associated with a severe disease course 
(6). STAT1 -7- mice exhibited a similar sustained, yet ineffec¬ 
tive, immune response. However, these clinical cases exhibited 
gene expression more characteristic of the type I interferon 
response and not of the T H 2, M2, and T H 17 cell expression 
observed in STAT1 -7- mice. These differences may reflect the 
choice of tissue used for microarray analysis (Cameron and 
coworkers used blood from patients whereas our analysis used 
lung tissue), but they could also be the result of the profound 
dysregulation associated with the loss of STATl-dependent 
transcription. However, similar T H 2 bias and macrophage dys¬ 
regulation characteristics have been observed in aged human 
populations (26, 59), members of which are more susceptible 
to pulmonary fibrosis (44) and exhibited greater mortality rates 
during the SARS outbreak (9). Based on these common fac¬ 
tors, STAT1 -7- mice may be a good model for investigation of 
the SARS-dependent initiation of fibrosis observed in aged 
populations. 

In conclusion, rMA15 infection of STAT1 -7- mice resulted 
in an aberrant T H 2 response, sustained viral infection, inflam¬ 
mation, and a profibrotic phenotype. This infection system 


serves as a model for the early and intermediate stages of 
ARDS and should be of value for the testing of therapeutic 
interventions to alleviate virus-induced ARDS. In addition, 
this model may be useful for studying the role of eosinophils 
and T h 17 cells in mediating persistent inflammation in re¬ 
sponse to viral infection. 
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